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ABSTRACT 

In this paper, a new efficient computational algorithm is presented for solving cyclic heptadiagonal linear 
systems based on using of heptadiagonal linear solver and Sherman-Morrison-Woodbury formula. The 
implementation of the algorithm using computer algebra systems (CAS) such as MAPLE and MATLAB is 
straightforward. Numerical example is presented for the sake of illustration. 
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1. INTRODUCTION 

The cyclic heptadiagonal linear systems often occur in several fields such as numerical 
solution of ordinary and partial differential equations, interpolation problems, boundary value 
problems, etc. [1,2]. In many of these areas, it is necessary to obtain the solution of the cyclic 
heptadiagonal linear systems. 

In this paper, we consider general cyclic heptadiagonal linear systems of the form 

HX =R, (1.1) 

where 
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X =(x 1 ,x 2 ,...,x n f, R =(R 1 ,R 2 ,...,R n f ,n >8. 

In paper [3], Karawia give an efficient symbolic algorithm to obtain the inverse of 
heptadiagonal matrix of the form (1.2) and then the solution of the systems (1.1) based on LU 
decomposition. There are many special cases of cyclic heptadiagonal linear systems (See [4- 

11]) 



Recently in [12], a new efficient computational algorithm is presented for solving nearly 
penta-diagonal linear systems based on the use of any penta-diagonal linear solver. In this 
paper we are going to compute the solution of a general cyclic heptadiagonal system of the 
form (1.1) without imposing any restrictive conditions on the elements of the matrix H in 
(1.2). Our approach is mainly based on getting the natural generalization of the algorithm 
presented in [12]. The development of a symbolic algorithm is considered in order to remove 
all cases where the numerical algorithm fails. 

The paper is organized as follows. In Section 2, symbolic computational algorithm for the 
solution of heptadiagonal linear system, that will not break, is constructed. In Section 3, 
Sherman-Morrison- Woodbury formula is given. An illustrative example is given in section 
4. In Section 5, Conclusions of the work are presented. 

2. Heptadiagonal linear solver 

In this section we shall focus on the construction of new symbolic computational 
algorithms for computing the solution of general heptadiagonal linear system of the form: 



H h X h ~ R h> 



(2.1) 



where 
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(2.2) 



Firstly we begin with computing the LU factorization of the matrix H h . It is as in the 
following: 



H h =LU 



(2.3) 



where 
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(2.4) 



(2.5) 



The elements in the matrices L and U in (2.4) and (2.5) satisfy: 



a. =\ 



d x if /' =1 

d 2 -f 2 g l Hi =2 

d^-e i z x -f i g 2 Hi =3 

d. -^ J -C i _ z -e i z i _ 2 -f i g i _ 1 ifi =4,5,...,m. 
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if i = 4,5, ...,m, 

if i = 1 
if i = 2 

if z = 3,4,..., w —1, 
if i =1 

if z = 2,3, w —2. 



We also have: 

dettf, =[]a, 



i=i 



Then, the solution is given by 
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At this point it is convenient to formulate our first result. It is a symbolic algorithm for 
computing the solution of a heptadiagonal linear system of the form (2.1). 

Algorithm 2.1. To compute the solution of a heptadiagonal linear system of the form (2.1), we may 
proceed as follows: 

Step 1: Set oii=di. If 0^=0 then a x = t end if. Set g^a^ z^A^ ki=A n _i/ a y f 2 =b 2 / a 1( e 3 =B 3 / a y a 2 =d 2 -f 2 
*gi. If a 2 =0 then a 2 = t end if. Set g 2 =a 2 -f 2 *Z!, f 3 =(b 3 -e 3 *g 1 )/a 2 , a 3 =d 3 -e 3 *z r f 3 *g 2 . If a 3 =0 then 
a 3 = t end if. 
Step 2: Compute and simplify: 
For i from 4 to m do 

e i =(B i -D i * gi- 3 /ai_ 3 )/ a,. 2 
fi=(b r Di* Zi_ 3 / ai_ 3 -ei*gi_ 2 )/ a hl 
Zi- 2 =Ai_ 2 -fi_ 2 *Cj_ 3 
gi-i=an-f|.i*z,.2-en*C|.3 



oti=(di- Di* C,J a|.3-e,z,.2-f,*gn) 
If a,=0 then c<j= t end if 
End do 

f n \ 



Step 3: Compute det H h = 



V 1=1 A=o 



Step 4: If det H h *0then do 

Step 5: Set Q hl =R hl , Q h 2=Rh2-f 2 *Qhi, Qh3=Rh3-e 3 *Qhi-f3*Qh2, 
Compute and simplify: 
For i from 4 to m do 

QhrRhi-D^Qhi.3/ a, 3 - e^Qh^-f^Qh,! 
End do 
Step 6: Set x hm =Q hm / 

a rrv x hm-l-(Qhm-l"gm-l* x hm)/ a m-l» x hm-2-(Qhm-2~gm-2 * x hm-l~ z m-2* x hm)/ °tm-2 i 

Compute and simplify: 
For i from m-3 by -1 to 1 do 

x hi = (Qhi"gi * x hi+l-Zi*Xhj + 2-Cj*Xhi + 3)/ . 

End do 

Step 8: Compute and simplify the solution: 
For i from 1 to m do 

Xhi = ( x hi)t=0. 

End do 

Else 

OUTPUT("The matrix H h is singular"); Stop. 

End If 

The new algorithm 2.1 is very useful to check the nonsingularity of the matrix H h . 



3. Sherman-Morrison-Woodbury formula[13] 

In this section, we are going to formulate a new computational algorithm for solving cyclic 
heptadiagonal linear systems of the form (1.1) based on the previous heptadiagonal linear solver. 
The heptadiagonal linear system of the form (1.1) can be written in the form: 



where 
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Thus (3.1) is equivalent to 



D n 
A n 

x n -2f ,x" = (x n _ l ,xj , R' = (R l ,R 2 ,...,R n _ 2 f ,andRHR„-M T - 



M x x'+V x" = R' 
U T x' + M 2 x" = R" 

Assume that M 2 is nonsingular. After elimination of x " from (3.2), we get the linear systems 



(3.2) 



M x' = R 



(3.3) 



where M =M, -VM^U 1 ,R=R' -VM 2 R ' 
If we applying the Sherman-Morrison-Woodbury formula to M, we obtain 

M' 1 =M 1 _1 +M;V (M 2 -U t M;V y x U T M; { and jc ' =M~ l R = y +m;V (m 2 -u t m;V y x u T y 

where y is the solution of M t y = R. It is clear that the solution x " can be found from the above 
formula by successive calculation of the expressions 

y = M 1 ~ 1 R, M;V , U T M;V , (M 2 -U t M~V ) ',and (M 2 -U T M[V Y l U T y . 
The main part of above calculations is in finding the first two expressions, which is equivalent to 
solving three (n _ 2)-by-(n _ 2) heptadiagonal linear systems with the same coefficient matrix IV^ and 
different right-hand sides. After finding ofx', we can get x" from the second equation of (3.2) by 
formula 

x"=M-\R n -U T x') 

At this point it is convenient to formulate our second result. It is a symbolic algorithm for 
computing the solution of a cyclic heptadiagonal linear system of the form (1.1) and can be 
considered as natural generalization of the symbolic algorithm 1 in [12]. 

Algorithm 3.1. To compute the solution of a cyclic heptadiagonal linear system of the form (1.1), we 
may proceed as follows: 

Step 1 : Find Mi , M 2 , U T , V, R',R \ and R = R ' -VM 2 R ". 

Step 2: Solve M x y =R,M l q l =v v mdM l q 2 =v 2 by algorithm 2.1, then obtain y and 

M-'v =(q l ,q 2 ). 

Step 3: Compute x' = y + (q v q 2 )(M 2 -U T {q x ,q 2 )y l U T y ,x" =M 2 l (R" -U T x'). 
Step 4: Compute the solution x = 



Three systems M x y -R,M x q x =v l5 mdM l q 2 -v 2 in algorithm 3.1 can be solved in parallel. 



4. An illustrative example 



In this section, we give simple numerical experiment to illustrate the effectiveness of our 
algorithms. 

Example 4.1. Consider the 10-by-10 cyclic heptadiagonal systems coming from [3] 
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Solution: The application of Algorithm 3.1 gives 
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Step2: y = 



2 2 2 
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M[V =(q x ,q 2 ) 



Step3: x ' = y + (q l ,q 2 )(M 2 -U T (q x ,q 2 )y x U T y =(1,2, 3, 4, 5, 6, 7, 8), 
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x"=M- l (R"-U T x') = (9, 10). 



Step 4: The solution of (4.1) is 



x = 



\ x J 



= (l,2,3,4,5,6,7,8,9,10) r . 



5. Conclusions 



In this paper, we derived a computational algorithm for solving the cyclic heptadiagonal linear 
systems. Since the algorithm uses the symbolic algorithm 2.1 method for the LU factorization of the 
heptadiagonal matrix, the factorization never suffers from breakdown, and this leads to the fast and 
reliable solution of cyclic heptadiagonal linear systems. The realization of the method needs O(n) 
operations. The algorithms are natural generalizations of some algorithms in current use. 
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